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Abstract We investigate the surface of the chiral phase transition in the three- 
dimensional parameter space of temperature, baryon chemical potential and mag- 
netic field in two different approaches, the field-theoretical Nambu-Jona-Lasinio 
i <~| (NJL) model and the holographic Sakai-Sugimoto model. The latter is a top-down 

Q-| approach to a gravity dual of QCD with an asymptotically large number of colors 

and becomes, in a certain limit, dual to an NJL-like model. Our main observation is 
D , that, at nonzero chemical potential, a magnetic field can restore chiral symmetry, in 

1 apparent contrast to the phenomenon of magnetic catalysis. This "inverse magnetic 

catalysis" occurs in the Sakai-Sugimoto model and, for sufficiently large coupling, 
in the NJL model and is related to the physics of the lowest Landau level. While in 
most parts our discussion is a pedagogical review of previously published results, 
we include new analytical results for the NJL approach and a thorough comparison 
of inverse magnetic catalysis in the two approaches. 



1 Introduction 

Two of the most important laboratories for the theory of strong interactions ex- 
hibit large magnetic fields: firstly, in non-central relativistic heavy-ion collisions the 
magnetic field perpendicular to the collision plane can be as high as 10 18 G f l |, and, 
secondly, in certain compact stars called magnetars the surface magnetic field is of 
the order of 10 15 G |2), while the application of the virial theorem suggests that in 
the interior the magnitude of the magnetic field might reach 10 18 G (3). Since this 
is comparable to the QCD scale Aqcd — 200 MeV [in natural Heaviside-Lorentz 
units, 10 18 G ~ (140 MeV) 2 ], the magnetic field in these laboratories might have a 
significant influence on properties governed by the strong interaction. For example, 
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in the case of heavy-ion collisions, the magnetic field might be responsible for an 
observed charge separation, which has been attributed to the so-called chiral mag- 
netic effect 13][5j[6)- On the other hand, in compact stars the structural composition 
of the star's interior, i.e., the equation of state, and transport properties could be 
affected . 

Also from a theoretical point of view these two physical systems present a great 
challenge. Both of them cover regions of the QCD phase diagram that are very 
difficult to study from first principles, since the large coupling strength prevents 
the application of perturbative methods. Relativistic heavy-ion collisions explore 
the phase diagram at low chemical potential and intermediate temperature (of the 
order of the QCD scale). This region is best tackled by lattice QCD, which has 
in the recent years been able to quantify the equilibrium properties in this regime 
H). F° r transport phenomena, however, lattice simulations are not well suited. 
Here, the application of the AdS/CFT correspondence [9|, a method developed in 
string theory, has contributed the celebrated result for the ratio of shear viscosity 77 
over entropy density s iflOl . The result, rj/s = 1 /An, is currently unrivalled by other 
methods and appears to agree very well with experimental data. Furthermore, this 
value has been conjectured to be a lower bound for all isotropic fluids 111 II . (This 
bound has been lowered in higher-derivative gravity duals 1 12), while in anisotropic 
fluids there appears to be no lower bound 1131 .") Unfortunately, the region of the 
phase diagram relevant for compact stars, where the temperature is low and the 
quark chemical potential is of the order of the QCD scale, is inaccessible for lattice 
simulations due to the so-called sign problem. Here one has to rely on extrapolations 
(down in density) from perturbative calculations or extrapolations (up in density) 
from nuclear physics or on suitable models, two of which will be of relevance for 
this article: the Nambu-Jona-Lasinio (NJL) model IT4l[T5l and the Sakai-Sugimoto 
model US El. 

In this review we focus on the effect of a background magnetic field on the 
chiral phase transition of QCD. The Lagrangian of QCD with Nf flavors exhibits 
an approximate global U(NAl x U(Nf)R symmetry group, which is broken down 
to a global SU(Nf)L x SU(Nf)R x U(1)l+r by the axial anomaly. At small tem- 
peratures and chemical potentials, i.e., in the hadronic phase, the chiral symmetry 
SU {Nf)L x SU(Nf)R is spontaneously broken to SU (Nf)L+R through the formation 
of a quark-anti-quark condensate. In this scheme the light mesons are understood 
as the (pseudo-)Goldstone modes corresponding to the broken generators of the 
symmetry group. By turning on a chemical potential one introduces an asymme- 
try between quarks and anti-quarks and thus exerts a stress on their pairing. As a 
consequence, one expects to eventually restore chiral symmetry^ 



1 At asymptotically large chemical potentials it is known from first principles that chiral symmetry 
is also broken, however via a different mechanism, namely by the formation of a diquark conden- 
sate in the color-flavor-locked phase |18||19| . Whether the hadronic phase is superseded by normal 
quark matter or by CFL or by some other color-superconducting phase is a matter of debate. Here 
we shall ignore color superconductivity. For the inclusion of color superconductivity in the context 
of the chiral phase transition in a magnetic field see Ref. 1201 . 
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The two models under consideration realize the implementation and breaking of 
chiral symmetry quite differently. In its original formulation, the NJL model was 
supposed to explain the mass of nucleons via chiral symmetry breaking. With the 
advent of QCD it was reformulated as a model of quarks ||2T1 l22l . It is a non- 
renormalizable field theory since it approximates the interaction of quarks by a 
four-point fermion interaction, and therefore the results of the model depend on 
the regularization scheme and on the UV cut-off that is used. Furthermore, the NJL 
model in its standard form lacks confinement. In the chiral limit, the Lagrangian of 
the NJL model is invariant under the same symmetry group as the QCD Lagrangian 
with massless quarks - the global SU(N C ) x U(N/)l x U(N/)r, where N c denotes 
the number of colors. The chiral symmetry is broken explicitly by a bare quark mass, 
which has to be sufficiently small compared to the momentum cut-off. Spontaneous 
breaking of chiral symmetry is then realized very similarly to the BCS theory of 
superconductivity ll23ll . which actually served as a guideline in the development of 
the model. Nevertheless, there are important differences between the condensation 
of Cooper pairs in a superconductor and the condensation of fermion-anti-fermion 
pairs. For instance, the presence of a Fermi surface in the former case implies the 
instability with respect to Cooper pairing for arbitrarily weak attractive interactions, 
while, as we shall see in the NJL model, there is a finite critical coupling that is 
needed to form a chiral condensate. We shall also see that the analogy becomes 
better for chiral condensation in a magnetic field. 

Our second model, the Sakai-Sugimoto model, is a top-down string-theoretical 
approach to a holographic dual of large-A^- QCD. It exploits a non-supersymmetric 
variation of the original gauge-gravity duality conjectured in J9) known as the Wit- 
ten model |24|. Sakai and Sugimoto introduced fundamental quarks in the chiral 
limit by placing a stack of Nf probe D-branes for the left-handed and anti-D-branes 
for the right-handed sector into the supergravity limit of the Witten model. Ac- 
cording to the holographic dictionary, the local gauge symmetry on these "flavor 
branes" translates into a global symmetry on the boundary, i.e., into the chiral sym- 
metry of the dual field theory. The question of whether the symmetry is intact or 
broken amounts to asking whether one can perform gauge transformations on D- 
branes and anti-D-branes independently or not, i.e., whether the D-branes connect 
with the anti-D-branes in the bulk. Thus the symmetry breaking mechanism in the 
Sakai-Sugimoto model is of geometrical nature. 

In order to understand what effects a magnetic field might have on the formation 
of the chiral condensate, let us recapitulate the general discussion found in ||25l . 
Calculating the chiral condensate in field theory amounts to calculating a fermion 
loop. Let the bare fermion mass be finite for the moment and regularize the UV 
divergence via a cut-off in some suitable scheme, e.g., Schwinger's proper time 
method. In the presence of a magnetic field one has to take Landau quantization 
of the transverse momentum of the charged fermions into account. It turns out that 
if one performs the chiral limit on the result, an IR singularity appears, which can 
be shown to originate from the lowest Landau level. As a consequence, a mass gap 
is dynamically generated in order to avoid this IR singularity. The precise form 
of the gap is of course dictated by the form of the interactions in the model under 
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consideration. This effect - termed magnetic catalysis - was first found in the Gross- 
Neveu model (26] |27l and later on in several NJL model calculations |25l |28l |29l 
[301 and in QED (31] as well as in holographic approaches ll32l l33l [34l [35l [36l 
[371 [38) ■ It also plays an important role in the context of graphene fl39l l40l . For 
QCD, it was found in a lattice calculation (however, with unphysical quark masses) 
that the critical temperature increases with the magnetic field fiTI . in accordance 
with magnetic catalysis. However, recently the Budapest- Wuppertal collaboration 
found (with physical quark masses) that the maximum of the quark susceptibility 
drops significantly at temperatures about 140 MeV under the influence of a magnetic 
field 142L i.e., the opposite effect was observed. It remains an open and interesting 
question what prevents magnetic catalysis to persist for larger temperatures in QCD. 

This article is mostly, but not exclusively, a review of existing work in the field- 
theoretical NJL and the holographic Sakai-Sugimoto model. In Sec. [2] we review 
the effect of a magnetic field on chiral symmetry breaking in the NJL model in a 
pedagogical way, starting from the simplest case without magnetic field. This sec- 
tion also contains several new, so far unpublished, aspects, for instance the analytic 
approximations and related discussions regarding inverse magnetic catalysis in Sec. 
12.2.31 After a pedagogical introduction to the Sakai-Sugimoto model in Sec. 13. II we 
discuss its limit of a small asymptotic separation of the flavor branes and map out 
the critical surface of chiral symmetry breaking in the parameter space of tempera- 
ture, chemical potential and magnetic field. This part is mostly a review of our own 
works l43l and fi4ll . with emphasis on the analytic approximations of the results 
and their comparison to the field-theoretical analogues. 



2 Chiral phase transition in the Nambu-Jona-Lasinio model 

We start from the standard Lagrangian of the NJL model (for an overview over the 
various NJL-type models see B31 ). 



We restrict ourselves to Nf — 1 ; pL denotes the quark chemical potential, m the bare 
quark mass, = + iqA^ the covariant derivative with q the electric charge of the 
quarks and the electromagnetic gauge potential = (0, —A). We shall work with 
imaginary (Euclidean) time x = —ix° compactified on a circle, the circumference of 
which is identified with the inverse of the temperature T. 

As a next step, we assume the pseudoscalar condensate to vanish, (xj/Ys W) = 0, 
and apply the mean-field approximation, 



if = y(iY fl D l 



m + jj.Yo)\j/ + G (xj/vf + {YY5Y) 2 



(1) 



(y/y/) 2 ~ -(v/v/) 2 +2(v/v/)v/v/. 



(2) 



We assume the quark-anti-quark condensate (y v) to be homogeneous and isotropic. 
For more general ansatze see |[46]|47], where a dual chiral density wave (a.k.a. chi- 
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ral spiral) has been discussed, and [481 for more general inhomogeneous phases. We 
define the constituent quark mass as 

M = m-2G(\j7y). (3) 

In the following we assume stationarity and apply the temporal gauge fixing con- 
dition. Therefore, the temporal dependence of the eigenfunctions of the Dirac op- 
erator is simply an exponential, e lfi) " T , with the fermionic Matsubara frequencies 
CO„ = (2n+ \ )nT. Then, the thermodynamic potential becomes 

= _T = (M-rnf _ T -jq, + g - g 
V 4G V T 

where Z is the partition function. The trace includes summation over the fermionic 
Matsubara frequencies and some suitable spectral decomposition e of the Dirac 
Hamiltonian, 

#D = rV(-iV -qX) + fM. (5) 

In analogy to BCS theory, the equation for minimizing the effective potential with 
respect to M is called gap equation, which reads 

t y° 

CWw) = Tr . (6) 

In the context of a background magnetic field, we shall also discuss the axial current. 
Its expectation value is given by 



<^> = 4Tr-^p^. 
V iat n + ii-e 



2.1 Chiral symmetry breaking without external fields 

Without external fields, the normalized eigenfunctions of the Dirac Hamiltonian are 
given (in a Weyl basis) by the momentum eigenfunctions 

y/ e , s , k (x) = -J=e ik - X ^= Vft - esk,^ s e^e k + eskY , (8) 

where £, s are two- vectors defined by the eigenvalue equation k • <7^ v = s% s with the 
usual Pauli matrices a', and where ee^ = ey/k 2 +M 2 with e = ± is the eigenvalue of 
the Dirac Hamiltonian, which in turn is two-fold degenerate with respect to s = ±. 
For the diagonal matrix elements of the gamma matrices in this basis we find 



&,* = «—. (rVk-y) =s. 

tfc \ / e,s,k 



(9) 
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From the second relation and Eq. (0 we conclude that the axial current vanishes, 
(A) = 0- The reason is the spin degeneracy in each state. This will no longer be 
true in the presence of a magnetic field, as we shall discuss in Sec. 12.2.31 We can 
compute the thermodynamic potential and, by inserting the first relation into Eq. 
©, the gap equation in the thermodynamic limit (at vanishing magnetic field B), 



&b=o = 



(M - to) 2 
4G 



2E 

e=± 



d 3 £ 



Tin 1 +e 



T 



(10) 



M — to =4G 



d 3 k M 



[l-f(e k -fi)-f(e k + ii)} , 



(11) 



where f(x) = \/{e x l T + 1) is the Fermi-Dirac distribution function. The (vacuum 
parts of the) momentum integrals are UV divergent and have to be regularized. Since 
the NJL model is non-renormalizable, all results, e.g., the magnitude of the gap 
and the order of phase transitions, will depend on the regulator as well as on the 
regularization scheme. We use the proper time regularization scheme ll49l . In this 
procedure, the integrand of divergent expressions is recast into so-called proper time 
integrals, 

7\ — a 1 f°° , „_i ^i^j-fJl 



(* 2 



r a) 



dTT"-'e- T ^+ fo ). 



(12) 



and one then performs the momentum integral before the proper time integral. The 
UV divergence of the momentum integral reappears at the lower bound of the proper 
time integral, which therefore has to be regularized. We set the lower bound to 1/A 2 . 
This yields the thermodynamic potential at zero temperature 



, 2A 2 (M-to) 



-A 2 (A 2 -M 2 ) e- M2 /A 2 +M 4 r (o, ^ 



2B(ii-M) 



tf(2, 2 -5M 2 ) + MUntt±^ 



(13) 



where r(a,x) is the incomplete gamma function, and the gap equation 



M — to = Mg 



-M 2 /A 



2 M 1 



M- 



A 2 " F % 2 



2Mg9(ii-M) 



A 2 



M 2 n+k F 
— rln 

A 2 M 



(14) 



where we have defined the the Fermi momentum kf = y/ [i 2 —M 2 and the dimen- 
sionless coupling constant 

GA 2 

* S 2^- (15) 



Inverse magnetic catalysis 



7 



For simplicity we shall discuss the chiral limit m = in the rest of the paper. In this 
case, M = is always a solution to the gap equation. 

For /i = 0, the gap equation further simplifies since the term <x 6 (ji — M) does 
not contribute. Then, after dividing Eq. ( [Pfl ) by M and g, its right-hand side is always 
smaller than 1 . Therefore, a nontrivial solution for M only exists if the dimensionless 
coupling constant g is larger than 1 . When it exists, this solution is preferred over 
the trivial solution, as one can verify with the help of the thermodynamic potential 

<m. 

In Fig. [T]we show the numerical solution for the gap equation as a function of jU 
for three different coupling constants larger than 1 (i.e., they all admit a nontrivial 
solution for fi = 0). For all couplings, there is a certain critical jj. where M goes to 



< 
S 



: T=0 


#=1.69 


g=go 




g=1.09 




1 




.1 





Fig. 1 The zero-temperature solu- 
tion to the gap equation for three 
different values of the coupling g. 
The thin dotted line is the line fl = 
M. The solution becomes multi- 
valued in the region /I > M for all 
couplings larger than go with go 
given in Eq. )19> . 



fi/A 



zero. By first dividing the gap equation by M and then setting M = 0, it is easy to 
show that this critical fi is given by 




(16) 



If and only if the solution is single-valued, this is the critical at which the (then 
second-order) phase transition to the chirally restored phase occurs. 

Above a certain coupling, the solution becomes multi-valued. The coupling 
where this qualitative change occurs can be computed as follows. By differentiating 
the gap equation with respect to fi we find 



dM 



2k F 



M 



(17) 



In accordance to the numerical plot, this derivative is infinite for M = 0. For all 
couplings for which the solution is multi-valued, there is another point where the 
derivative is infinite, which is given by the second pole of the denominator, 
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r(0,M 2 /A 2 ) 

H =Mcosh 2 -• (18) 

We can now ask for the value of g at which this point coincides with fio(g) for 
M —> 0. The resulting equation then yields the coupling where the multi-valuedness 
sets in. We find 

80 = -^-^ 1.390, (19) 

2 

where Je is the Euler-Mascheroni constant. In the regime 1 <g < go the chiral phase 
transition is second order and takes place at Ho(g). 

For couplings larger than go the transition is first order and has to be determined 
numerically. It turns out that the branch with a positive slope is always energetically 
disfavored. Therefore, in terms of Fig. Q] the preferred solution follows the horizon- 
tal line M(jU = 0) and, for all multi-valued cases, jumps to zero at a certain chemical 
potential. Whether (and how far) the preferred solution follows the curve into the re- 
gion jj. > M depends on the coupling. We find numerically that for couplings below 
(above) g ~ 2.106 it does (doesn't). This is a first example of the nontrivial effect 
of jJ, on the preferred phase: it is not always the phase with the largest dynamical 
mass that is favored. In more physical terms, for couplings above g ~ 2.106 the chi- 
rally broken phase with vanishing quark density is directly superseded by the quark 
matter phase, while for smaller couplings there is a region of finite density between 
these two phases. Since for g > 2. 106 there are no complicated effects of the quark 
density, we can write down a very simple expression for the free energy difference 
between the broken phase and the restored phase, evaluated at the solution of the 
gap equation (and using M <C A), 

All = -$*(l-l) + J* t (20) 



167T 2 V gj I2n 2 

with Mq being the (non-analytical) solution to the gap equation for jj. = 0. This 
result is very intuitive: the first, negative, term is the condensation energy, i.e., the 
energy gain from the chiral condensate, while the second, positive, term corresponds 
to the energy costs for pairing which must be paid because the chemical potential 
has separated fermions from anti-fermions. When the costs exceed the gain, chiral 
symmetry is restored. This determines the phase transition line. Below we shall 
derive the analogue of this strong-coupling free energy difference in the presence 
of a background magnetic field, see Eq. d40b . We summarize our discussion of the 
chiral phase transition at B = T = in Fig. [2] 

For nonzero temperatures, we need to solve the gap equation ( fTTT ) [with the regu- 
larization of the vacuum part shown in Eq. (TBI) ! numerically. The result for various 
temperatures and a large coupling (larger than that of point c in Fig. O is shown 
in the left panel of Fig. [3] In general, the temperature decreases the gap. Moreover, 
the temperature can also change the order of the chiral phase transition by removing 
the multi-valuedness of the solution to the gap equation. The critical temperature of 
the chiral phase transition in the T-jj, phase diagram is shown in the right panel of 




Fig. 2 The phase diagram at T = in the /l-g-plane. Dashed lines indicate second-order, solid lines 
first-order phase transitions. In the shaded region chiral symmetry is restored OfS). The points a, 
b and c correspond to (fl/A,g) = (0, 1), (e _1fe / 2 /2,go)> an d (0.542,2.106), respectively, with go 
given in Eq. )19t . Between points a and b the transition line is given by (UoG?) from Eq. U6\ . The 
dashed line between a and c indicates the onset of a finite quark number density n q within the 
chirally broken phase (^Sb). 




Fig. 3 Finite-temperature effects on the chiral phase transition in the NJL model. Left panel: the 
gap as a function of the chemical potential for a given coupling strength and different values of 
temperature. Right panel: the phase diagram in the ^-T-plane for the same coupling. The (dashed) 
second-order phase transition line is given by the analytic expression i2l\ . 



Fig. [3] The critical point moves towards higher temperatures with increasing cou- 
pling. If the phase transition is second order it is possible to find a closed form for 
the critical temperature. To this end, one divides Eq. (fTTT i (with m = 0) by M and 
then sets M = in the remaining equation. Then, solving for T yields the critical 
temperature 

A 




(21) 
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2.2 Chiral symmetry breaking in the presence of a magnetic field 



2.2.1 Structure of the fermion states in a background magnetic field 

Let us consider a homogeneous background magnetic field B = (0, 0, B) by choosing 
the Landau gauge fixing condition with A = (— yB, 0,0). Within this ansatz, the 
eigenfunctions of the Hamiltonian are proportional to exp [i ( co„ T + k x x + k z z)] . Using 
this, we split the Dirac Hamiltonian in a longitudinal and a transverse part with 
respect to the direction of the magnetic field, Hp =H^ + Hj, where 

1 0\ [0 a f 



Hl = ffh + fM, H T = sgn(^) ( q 1 ) ® ( a ' ' (22) 



with 



^ +s8 " ta) '^b<-^ <23) 

We see that a is the annihilation (creation for q < 0) operator of the quantum me- 
chanical oscillator, which gives rise to the Landau quantization of the energy spec- 
trum of a charged fermion moving in a background magnetic field. For q > 0, the 
orthogonalized eigenfunctions of the Hamiltonian are given by 



& i{k z z-\-k x x) j 



y/K k _J - sk zy JS k _J ~ SK k _t (4 |^ - 1) 

es i 



^Kk-j + sk zy /e k .j + sK k _j {£, \l) 
\ey/K k .J - *k zy /e kz j + SK k _J (%\t-l)J 

where I = 0, 1,2,3, . . . denotes the Landau level, where 

1/4 



(24) 



«w=^«(ir)^ H4|sp/2H '(^)' < 25 » 



1) =0, and 



£ M = ^kJ+M^2^\Bi, K kz/ = ^kj+2\q]m. (26) 

Here, Hi is the £th Hermite polynomial and L, the length of a box with volume V in 
the ;th direction. In order to obtain the eigenfunctions for the case q < 0, one simply 
replaces \i) with (5, \£ — 1) and vice versa. For the diagonal matrix elements of y° 
and y°y 3 75 we find 

t,,., = e^- , (rVrs W - sgn(9)^- ■ (27) 

%i KjCr,t 
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From Eq. (l24t we see that in the lowest Landau level (LLL) I = only the sgn(g) s = 
1 -states survive, which are also eigenstates of the spin operator^ = y°y 3 75 as well 
as zero-eigenmodes of Hj. Therefore, the dynamics of the LLL becomes effectively 
1 + 1 -dimensional. Moreover, in the limit M — >• for sgn(q)e k z > (< 0) these 
states are right- (left-) handed only. This is an indication that the magnetic field 
induces an axial current l50l . More precisely, due to the sum over s in the axial 
current (0, the relation ( l2Tb shows that only the LLL level contributes. Due to the 
sum over e there can only be a finite contribution if jj. 7^ 0. Since we have put 
the fermions into a box with volume V = L x L y L z , the range of y is restricted to 
[— L y /2,Ly/2] and therefore k xjnax — k xjn i n = L y \q\B since we have absorbed k x into 
the new coordinate Hence, because of Ak x = 2%/L x , each energy level for given 
e, kz, s and I has a degeneracy of L x L y \q\B/ (2k). In two cases the result for the 
axial current along the magnetic field can be given in closed form, 

M = 0,Vr: Ul) = ^, (28) 

T = 0,VM<»: UD = qB ^T^ . (29) 

The prefactor \q\B/(2n) found by phase space considerations has a very special 
role here. It is the difference of the number of zero-eigenmodes of Hj with s = 1 
and s = — 1 respectively. This is a topological result since it is given by the index of 
each 2x2 block of Hj, which in turn is linked to the Euclidean chiral anomaly in 
two dimensions via the index theorem. Furthermore, the first result is independent 
of T which is a special feature of massless 1 + 1 dimensional fermions and hence 
reflects the effective dimensional reduction. 



2.2.2 Magnetic Catalysis 

Let us return to chiral symmetry breaking, now in the presence of a magnetic field. 
The thermodynamic potential and the gap equation read 



\q\B 



4G 2k 



E I> 

e=±e=o 



dk z 
,2k 



Tin 1 +e~ 



(30) 



M =2G 



\q\B 
2k 



1=0 



dk, M 



I> / + , (3D 



where a^ = 2— 5qc. Comparing with the corresponding B = expressions in Eqs. 
( [Tol l and (fTTl i. we see that the effect of the magnetic field is to replace —> 1 and 
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Using again proper time regularization, the thermodynamic potential at vanishing 
temperature becomes 



£2t=q — A 



fi=T=B=0 - 



{qBf 

2n 2 



r* 4 



T (3-21nx) + f(lnA_ 



>(*) 



\q\B 



^ 6 (M -M) £ a e (vk F ,c - M) In 



F.I 



Mi 



(33) 



Here, Qh=t=b=o is the vacuum part from Eq. ( fT3l l. y/"' the n-th polygamma 
function (analytically continued to negative values of «), we have abbreviated 
x = M 2 /(2\q\B), and 



M e = JM 2 + 2\q\Be, k 



F,e : 



H 2 -M 2 
~2^\B~ 



(34) 



Different regularization schemes - compare for instance with flSTl , where dimen- 
sional regularization is used - only differ in the B — result and in (divergent) 
terms that depend on B but are constant in M, which are omitted. The latter can be 
viewed as renormalizing the energy content coming solely from the magnetic field. 
The corresponding gap equation is 



M = Mg 



A 2 



A 7 



-2Mg 



\q\B 
A 2 



1 \ . N 1 

— — x J mx +x + lnf" (x) — — In 2n 



\q\B ^ U+k F e , 

£=0 £ 



(35) 



Let us first consider the case ji = 0, i.e., we can ignore the terms oc 9(jJ. —M) in Eqs. 
( f33b and d35l ). For small coupling g <C 1, the dynamical mass squared will be much 
smaller than the magnetic field, M 2 <C \q\B. Then, with M< A, the gap equation 
becomes 

2\q\B, 



1 



A 2 



■In 



\q\B_ 
7tM 2 ' 



(36) 



Now, there is a nontrivial solution for arbitrarily small g. This is in contrast to the 
case B = where chiral symmetry can be broken only for g > 1 . The solution is 
obviously 



M : 



MB — r~r~ — 
— — — Q WO. 

K 



(37) 



This qualitative effect of the magnetic field on chiral symmetry breaking was termed 
"magnetic catalysis" in |28| and was since observed in numerous different models. 
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Fig. 4 Effects of magnetic catalysis on the dynamical mass M and the critical temperature. Left: the 
gap at T = = for different couplings. The lowest coupling shown corresponds to a subcritical 
coupling at B = 0, i.e., its nonzero value is solely induced by B. Its behavior at small B is given 
by the exponential in Eq. )37t . Right: the critical temperature for chiral symmetry restoration as a 
function of B. 



Interestingly, as already mentioned in the introduction, this effect stems mainly from 
the physics in the LLL. In order to show that, one omits all contributions from I > 
in (l3lT l and cuts off the momentum integral at yJ\q\B/4n, since below that cut- 
off the LLL dominates. Then, one obtains exactly the result 071 ). Furthermore, the 
logarithmic IR singularity in d36l l regulated by the dynamically generated mass is 
precisely due to the LLL contribution and its 1 + 1 dimensional nature. The form of 
the gap in the weak coupling limit is reminiscent of the BCS gap in a superconductor 
ll52ll . In both expressions for the respective gap the relevant density of states appears 
in the denominator of the exponent. Here it is the density of states of the massless 
fermions at £k z ,e=o = 0, whereas in the BCS gap it is the density of states at the 
Fermi surface. In both cases the dynamics is essentially 1 + 1 -dimensional. While 
in BCS theory this effective dimensional reduction is a consequence of the Fermi 
surface, here it is provided by the magnetic field. Note that the dimensional reduc- 
tion is not in conflict with the Mermin-Wagner-Coleman theorem that states that no 
spontaneous symmetry breaking can occur in 1 + 1 dimensions. The reason is that 
the Nambu-Goldstone modes are neutral, and hence their motion is not restricted 
by the magnetic field. At extremely large magnetic fields the internal structure of 
these modes can be resolved which might invalidate this argument ll53l . 

We show the numerical solution of the gap equation for various coupling strengths 
for T = ji = in the left panel of Fig. |4] Magnetic catalysis also manifests itself in 
the critical temperature for chiral symmetry restoration, which, at ji = 0, monoton- 
ically increases with increasing magnetic field, see right panel of Fig. [4] 

2.2.3 Inverse magnetic catalysis 

We now include the contributions from a nonvanishing chemical potential ji. First 
we discuss the case of weak coupling which corresponds to M 2 <C \q\B. Since the 
chiral phase transition can be expected to occur at chemical potentials of the order 
of the mass gap, we may thus also assume ji 2 <C \q\B (we are not interested in 
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the physics far beyond the phase transition). As a consequence, we can employ the 
lowest Landau level approximation, i.e., drop the contribution of all higher Landau 
levels. Then, from Eq. (|33l l we conclude that the difference of the thermodynamical 
potentials of the chirally broken phase and the quark matter phase is 



Again, we find a very interesting analogy to superconductivity: the resulting ex- 
pression is exactly the same as for a BCS superconductor with mismatched Fermi 
momenta - first discussed by Clogston l54l and Chandrasekhar ll55l - after M is 
replaced with the superconducting gap A, \q\B with the average Fermi momentum 
(squared) of the constituents of a Cooper pair, and fi with the difference of the 
respective Fermi momenta. (Note that again the degeneracy factor of the LLL emu- 
lates the role of the Fermi surface.) 

To discuss the meaning of AQ for the chiral phase transition, let us first consider 
the case of a fixed magnetic field B and start from /i = 0, i.e., in the chirally broken 
phase. Upon increasing /I, we will reach the point jj, = Mj \/2 where AQ. changes 
its sign and thus the phase transition to the chirally restored phase occurs. This 
point is, in the context of superconductivity, called the Clogston limit. It occurs 
before the second term has a chance to contribute since still jj, < M. Now, more 
importantly for our purpose, let us again start in the chirally broken phase, i.e., from 
A £2 < 0, but now we increase the magnetic field at fixed ji (as we have just seen, 
for the discussion of the phase transition we may assume ji <M and thus ignore the 
term °= 9(ji —M)). Since we have started from a negative jx 2 — M 2 /2, increasing 
the magnetic field can only make AQ more negative because the dynamical mass 
increases with B. Consequently, the magnetic field only brings us "deeper" into the 
chirally broken phase. This is what we have expected from magnetic catalysis. 

However, as we will now explain, for g > 1 and finite chemical potential this ex- 
pectation is incorrect. We shall rather find that, for intermediate values of the mag- 
netic field, an increasing magnetic field does restore chiral symmetry. Let us, to this 
end, first discuss the numerical solution of the gap equation, see Fig. [5] Due to the 
sum over the Landau levels, the gap exhibits the well-known de Haas-van Alphen 
oscillations. Similar to the behavior found for B = 0, only the branches with a neg- 
ative slope of M(jj.) can be energetically preferred. Depending on the specific value 
of g there might be several phase transitions within the gapped phase into regions 
with jj. > M, i.e., with a finite quark number density, before entering the restored 
phase M = 0. In general it is also possible that the order of the phase transition into 
the restored phase oscillates between first and second order upon varying B: in the 
example shown in the plot, at vanishing magnetic field the phase transition is first 





via gap equation 



Inverse magnetic catalysis 



15 



0.6 



< 



; T=0, g=1.69 






.. lq|B/A*=0 




-j* 




/ \ f 




|q|B/A 2 =0.19; \ \ 




■ f\jq|B/A 2 =0.13 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

ft/A 

Fig. 5 The zero-temperature dynamical mass as a function of the chemical potential for different 
values of the magnetic field. For the lowest nonzero value of \q\B shown (solid line), Landau level 
oscillations can be seen. The magnetic field for the two other curves (dashed and dashed-dotted 
lines) is sufficiently large to suppress all Landau levels except for the lowest. 



order, while at \q\B/A 2 = 0.13 it is second order and at \q\B/A 2 = 0.19 again first 
order. We also see that the dashed (blue) curve for the lower magnetic field reaches 
farther in the ji direction than the dashed-dotted (black) curve for the larger mag- 
netic field. This is a surprise from the point of view of magnetic catalysis: it seems 
to indicate that the critical chemical potential for chiral symmetry breaking can de- 
crease with increasing magnetic field. We discuss this "inverse magnetic catalysis" 
in more detail now. 

To this end, let us consider the "cleaner" case of sufficiently large couplings 
where symmetry restoration happens in the region jj. < M for all magnetic fields. 
In this case, oscillations of the critical line in the phase diagram originate solely 
from the restored phase (not from the solution of the gap equation), and the phase 
transition is always first order. The numerically obtained phase diagram for such a 
case is shown in Fig. [6] From the arguments in the previous subsection, one might 
have expected that magnetic catalysis leads to a monotonically increasing critical 
chemical potential as a function of B (just like the critical temperature in the right 
panel of Fig. [4]). However, this is not the case: there is a region in the phase diagram 
where, upon increasing B at fixed /x, chiral symmetry is restored, in contrast to the 
weak-coupling case discussed below Eq. ( 1381 ). 

In order to understand this phenomenon, let us derive an analytic expression 
for A£2, analogous to the weak-coupling case. As discussed, for the given large 
coupling, the solution to the gap equation is simply given by the jj. =0 solution. For 
small magnetic fields, \q\B -C M , we can expand the solution up to second order in 
the magnetic field, 



M~M 



6M 4 r(0,M ( 2/A2) 



(39) 
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Fig. 6 Zero-temperature chiral phase transition in the plane of magnetic field and quark chemical 
potential at a rather large value of the coupling constant such that the phase transition is first order 
for all magnetic fields. (For smaller values the shape of the transition line is similar, but the order 
can vary between first and second.) Apart from oscillations at small B due to higher Landau levels 
in the chirally restored phase, the critical chemical potential decreases up to qB/A 1 ~ 0.5, see 
explanation in the text. The dashed-dotted line is the approximation to the phase transition line 
from Eq. <40t . 



where Mq is the solution for ji = B = 0. Inserting this solution into Eq. (1331) . we 
obtain the free energy for the chirally broken phase up to second order in B. The free 
energy for the chirally restored phase is, although we can set M = 0, complicated due 
to the sum over Landau levels. Let us therefore ignore the higher Landau levels. This 
seems to contradict our assumption of a small magnetic field which we have made 
for the chirally broken phase. Nevertheless, we shall see that the phase transition 
line obtained from this approximation reproduces the full numerical line in a region 
of intermediate magnetic fields. Since this is exactly the region where the "back 
bending" of the phase transition line is most pronounced, this serves our purpose 
to capture the main physics of the inverse magnetic catalysis. With Mq -^i A, the 
resulting free energy difference is 



AO. ~ -- 



Ml A 2 
16k 2 



1 -■ 



1 



\q\B. 



(qB) 



i r 



4k 2 ^ 24k 2 



l-12C'(-l)+lri 



2\q\B 



(40) 



(This is the generalization of Eq. (f20b to nonzero (but small) magnetic fields.) This 
expression allows for a qualitatively different phase transition line compared to the 
weak-coupling limit ( f38l > for the following reason. The term linear in B corresponds 
to the cost in free energy to form a fermion-anti-fermion condensate at finite ji. 
Importantly, this cost depends not only on ji, but also on the magnetic field. This 
is also true at weak coupling. However, in that case, the gain from the condensa- 
tion energy was also linear in B. This is different here: now, if we start from the 
chirally broken phase, i.e., from AQ < 0, increasing the magnetic field can lead to 
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a change of sign for AQ and thus restore chiral symmetry. This is what we have 
termed inverse magnetic catalysis in P31 . In this reference, we have also explained 
that the physical picture can be understood once again in analogy to superconductiv- 
ity, where, in the presence of a mismatch in Fermi momenta, it is useful to think of a 
fictitious state where both fermion species are filled up to a common Fermi momen- 
tum. Creating such a state costs free energy which may or may not be compensated 
by condensation. The point of inverse magnetic catalysis is that creating such a ficti- 
tious state (where fermions and anti-fermions are not separated by jj.) becomes more 
costly with increasing B, while B still enhances the dynamical gap due to magnetic 
catalysis. The magnetic field thus plays an ambivalent role by counteracting its own 
catalysing effect. 

This effect was first observed in the NJL model in |56l at T = and in 1571 
for the full three dimensional T-ji-B parameter space, and has been confirmed in 
various other calculations |20] |58] [59] |60] EU |62] |63] . Only for sufficiently strong 
magnetic fields the system enters the regime where magnetic catalysis is dominant. 
Typical fits of the model-parameters yield a cut-off of the order of a few hundred 
MeV B31 . Translating this into a scale for the magnetic field shows that the regime 
of magnetic catalysis is beyond the magnetic field strength expected in compact 
stars, and thus, if there is any observable effect of the magnetic field for the phase 
transition between hadronic and quark matter, it is inverse magnetic catalysis. 



3 Chiral phase transition in the Sakai-Sugimoto model 
3.1 Introducing the model 

The model discussed in this section is based on the conjecture that particular 
strongly coupled quantum gauge theories are equivalent to certain classical gravita- 
tional theories in higher dimensions. In the context of string theory, the first realiza- 
tion of this holographic principle known as AdS/CFT correspondence was proposed 
by Maldacena J9j- In a nutshell, it utilizes two different limits of describing so- 
called D-branes, which are dynamical objects in string theory that impose Dirichlet 
boundary conditions on the endpoints of open strings. On the one hand, a stack of 
N c D-branes hosts a maximally supersymmetric U(N C ) gauge theory coming from 
the massless excitations of open superstrings; on the other hand, the stack of D- 
branes is a massive object that curves space-time by coupling to gravitons - coming 
from the closed strings - with the strength X <*; g s N c , where g s denotes the string 
coupling. Now, let N c — > °° and keep X fixed. In the limit A < 1, gravity decouples 
from the open strings, whose low-energy effective theory is given by the mentioned 
U(N C ) super Yang-Mills theory. In the case of D3-branes, this gauge theory is four- 
dimensional. In the opposite limit, X 1, the stack of D-branes back-reacts strongly 
on the background. Gravity far in the asymptotic region also decouples from the sys- 
tem due to the gravitational red shift. Therefore, one can zoom in to the near-horizon 
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region of the space-time, which in the case of D3-branes is given by AdSs x S 5 . The 
idea behind the AdS/CFT duality is that the classical (super-)gravitational descrip- 
tion is fully equivalent to the quantum theory of the large-A^, large A limit of the 
super- Yang-Mills theory. This particular gauge/gravity duality, which has passed 
many nontrivial tests, has since been greatly generalized and also been used in the 
form of phenomenological (bottom-up) models. 

The Sakai-Sugimoto model lfl6l [171 is a string-theoretical top-down approach 
to large-A^- QCD. It is based on a proposal for a holographic dual of a non- 
supersymmetric large-A^- Yang-Mills theory in four effective dimensions by Witten 
ll24ll . In contrast to the original AdS/CFT correspondence, the background is pro- 
vided by the gravitational field of a stack of D4-branes. The dual field theory now is 
4+ 1 -dimensional since this is the dimension of the world volume of the D4-branes. 
The extra dimension is compactified on an S 1 and thus breaks supersymmetry on 
the field theory side: by imposing anti-periodic boundary conditions on the adjoint 
fermions, they obtain a mass of the order of the inverse radius of the S , called 
Kaluza-Klein mass Mkk- At one loop level, also the adjoint scalars become mas- 
sive. Hence, by choosing the radius of the extra dimension small enough and by 
restricting to low energies, one effectively breaks supersymmetry and effectively 
reduces the number of dimensions to 3 + 1 . However, there is a price to pay for in- 
troducing the extra dimension: in order to justify the supergravity approximation for 
the D4-brane background, the five-dimensional (dimensionful) 't Hooft coupling A5 
has to be large compared to Mj^£. This corresponds to a large four-dimensional (di- 
mensionless) 't Hooft coupling A = Xs/(2kM^). In this case, however, the mass 
gap of the field theory is of the same order as Mkk an d thus the Kaluza-Klein modes 
do not decouple. Only in the opposite limit A <C 1, where string corrections are im- 
portant and which thus is inaccessible, the Kaluza-Klein modes do decouple and the 
theory becomes dual to large-A^ c QCD in 3+1 dimensions (at small energies below 
the Kaluza-Klein scale). It has nevertheless turned out that the classical gravity limit 
of the D4-brane background is a remarkably useful tool for understanding certain 
nonperturbative properties of (large-A^-) QCD. 

An important property of the Sakai-Sugimoto model is the existence of a Haw- 
king-Page transition between a soft-wall and a black hole background, which en- 
codes a confinement-deconfinement transition. This feature can be understood either 
from power counting in N c of the corresponding thermodynamic potentials of the 
gravity backgrounds or by studying the dual to the Wilson line. Confined and de- 
confined phases correspond to two different geometric backgrounds which are, in 
coordinates made dimensionless by dividing by the curvature radius R, given by 
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IF ^ I 1 — t<KK i 

M = { " 3 , h d {u)={ M 3 , h c {u) = \ M 3 ' (42) 



deconf. 



and 



u ™ = {T)m = v ut= \t)w (43) 

Here, j3^ 4 is the period of JC4 - the coordinate of the additional S 1 - necessary to 
prevent a conical singularity atu = mkk in the confined phase. The curvature radius 
is related to the Yang-Mills coupling gYM by 

ZM K K 

where ^ 2 = a' is the squared string length. In the analytic continuation to Euclidean 
signature, time is also compactified to a circle with circumference f} t = T , anal- 
ogously to finite temperature field theory. Increasing the temperature shrinks the 
Euclidean time circle. At the point where the circumference of the time circle 
and the extra dimensional circle match, the Hawking-Page transition takes place. 
Apart from the metric field the Witten model also contains a nontrivial dilaton and 
Ramond-Ramond (RR) flux background given by 

e * = »3/ V F4 JMpk in „ ,45, 

where Q.^ is the volume of the 4-sphere. 

Sakai and Sugimoto introduced fundamental quarks by placing two stacks of Nf 
D8-branes with opposite orientation into the background in the so-called probe limit 
Nf <C N c , i.e., back-reactions on the geometry are neglected. In the asymptotic re- 
gion u — y °° the two stacks of D-branes are separated on the Kaluza-Klein circle. 
In the original model they reside at antipodal points. In the bulk, the D-branes are 
space filling in the field theory directions, x^, as well as in the S 4 and are specified 
by an embedding function in the u-xa, subspace. Before going to the gravity de- 
scription of the D4-branes one can interpret the underlying string picture as follows: 
strings connecting the D4 with the D8-branes carry one flavor and one color index, 
hence representing (massless) quarks in the fundamental representation, whereas 
strings stretching between D8-branes represent mesons. The local symmetry of the 
U(Nf) x U(Nf) gauge theory supported on the world volume of the stacks of D8- 
branes translates into a global symmetry via the holographic dictionary, which is 
interpreted as the chiral symmetry of the field theory. In the confined background, 
the two stacks of D8-branes are forced to join at m^k where the additional S l degen- 
erates and therefore form a single stack with gauge symmetry U(Nf). On the field 
theory side, this reflects the chiral symmetry breaking mechanism. One can use a 
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diagonal subgroup of the full symmetry group to introduce chemical potentials and 
electromagnetic quantities such as an external, non-dynamical magnetic field. Usu- 
ally the gauge is chosen such that for example the asymptotic value of the zeroth 
component of the Abelian gauge field is identified with the quark chemical poten- 
tial. Due to the probe limit, the deconfinement transition is not affected by a finite 
chemical potential, trivially leading to a phase diagram in the plane T-ji similar to 
the one discussed for large-A^ QCD in l64l . 

The low-energy effective theory describing the open string fluctuations is a non- 
Abelian Dirac-Born-Infeld (DBI) theory on the probe branes; calculating the fluc- 
tuations of the gauge field corresponds to calculating the meson spectrum. Indeed, 
after fitting the value of the 't Hooft coupling A and Mkk to the rho meson mass 
and the pion decay constant, the spectrum matches experimental data nicely. The 
mode expansion used in the calculation of the meson spectrum can also be used to 
link the Sakai-Sugimoto model to the Skyrme model. Apart from the DBI action, 
the dynamics of D8-branes in a background with nontrivial RR-flux is governed by 
a Chern-Simons (CS) action, since the D8-brane is magnetically charged under that 
flux. This contribution allows for introducing baryons in the model and is related to 
chiral solitons in the Skyrme model. Therefore, the full action reads 

S= Sdbi + Scs 

= T S f drd*x e~*Tty/\ det( gmn + 2na'& mn ) \ + ^ / C 3 Tr(27ra'^) 3 , (46) 

JDS O JDS 

where Ts is the D8-brane tension, and dCi = F4. Usually one integrates the last term 
by parts, omitting all boundary terms, to obtain a gauge-variant action where the 
RR-flux couples to the Chern-Simons 5-form. 

The Sakai-Sugimoto model also has a connection to the NJL model. In the "de- 
compactified" limit where the asymptotic coordinate distance between the D8- and 
anti-D8-branes is much smaller than the radius of the extra compactified dimension, 
L <C the Sakai-Sugimoto model is dual to a non-local NJL model 11651 . As 
a consequence, in the scenario with broken chiral symmetry, the D8-branes now in 
general join at uq > mkk- The difference uq — mkk is commonly interpreted as the 
constituent quark mass within a meson, which is realized as a string with both end 
points attached to the tip of the joined D8-branes hanging down to the bottom of the 
geometry. With a sufficiently small asymptotic separation of the flavor branes, it is 
also possible to find an energetically preferred phase with broken chiral symmetry 
in the deconfined background l66l . see Fig. [7] The resulting phase diagram at finite 
chemical potential was first discussed in l67l . By reducing L compared to M^, 
the temperature range where the system is confined becomes small compared to the 
temperature range governed by the deconfined and chirally broken phase. Eventu- 
ally, the resulting phase diagram resembles the NJL result (where no confined phase 
is present) shown in the right panel of Fig. [3] Consequently, the Sakai-Sugimoto 
model allows for interpolating between a non-local NJL model (L <C M^) and - 
modulo the above mentioned caveats - large-A^ t QCD (L = kM^). In the former 
limit, the flavor D8-branes do not probe deeply the background geometry produced 
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by the color D4-branes (which corresponds to neglecting gluon dynamics), while in 
the latter the gluons dominate. 





U = oo 
u = u 

11 = 11, 

Fig. 7 The chirally restored (left) and chirally broken (right) phases of the non-antipodal Sakai- 
Sugimoto model in the deconfined background. The calculation reviewed here determines which 
of the two D8-brane embeddings is favored as a function of temperature, chemical potential, and 
magnetic field for a small asymptotic separation L. Only in that limit (in which the dual field theory 
resembles the NJL model) does the chiral phase transition in the probe brane approximation depend 
on chemical potential and magnetic field. 



The effect of a homogeneous background magnetic field has first been considered 
in ll68ll . Shortly thereafter, the effect on the critical temperature for chiral symmetry 
restoration at vanishing chemical potential has been analyzed fl69l . Like in the NJL 
result from Fig. [4] the critical temperature increases with the magnetic field, which 
shows that the Sakai-Sugimoto model exhibits magnetic catalysis. Finite chemical 
potentials have been introduced together with a magnetic field in ITUl in the original 
Sakai-Sugimoto model. The deconfined, chirally symmetric phase was discussed in 
||7T1| . where a magnetic phase transition within the symmetric phase was found that 
is reminiscent of a transition to the lowest Landau level. The full phase diagram in 
the parameter space T-ji-B in the deconfined phase was presented in our work l43l . 
In particular, the effect of inverse magnetic catalysis effect was found and discussed 
in this reference. 

The Sakai-Sugimoto model can be developed further to include homogeneous 
baryonic matter, made from point-like approximations to the solitonic baryons men- 
tioned above l72l . Applications in the context of a background magnetic field have 
been studied in the confined l73l and deconfined (44] backgrounds, the latter study 
investigating the effect of baryonic matter on inverse magnetic catalysis. Here we 
shall mostly focus on the case without baryons, and only at the end of Sec. 13.61 
briefly review their effect on the phase diagram. 



3.2 Equations of motion and axial current 

In terms of an embedding function X((k) for the D8-branes, the induced metric reads 
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ck 2 / 1 \ 

-^P = u 3 ' 2 h d dt 2 + u i l 2 8 i jdoc i Ax i + m 3 / 2 ( x}h c + j-^J du 2 + u l/2 dQ% , (47) 

where the prime denotes derivative with respect to u. We work with one flavor, 
Nf = 1, and for the (dimensionless) 1/(1) gauge field we choose the ansatz 



2%a! 



a = 



1? 



dx M =flo(w)df + &r dbr+a3(«)cbr, (48) 



where b = Ina'B denotes the magnitude of the dimensionless magnetic field 
strength. Note that the necessity of introducing the third component of the gauge 
field, which is P-odd, is due to the coupling to ciq and b via the (P-odd) CS-action. 
We denote the asymptotic values of the gauge field by 

2na' I 
aoH = J u = jU? __ ; a 3 (oo)=j, x A (°o) = -, (49) 

where jx q is the dimensionful quark chemical potential, and £ = L/R is the dimen- 
sionless asymptotic separation of the flavor branes. The boundary value of 03 can be 
shown to correspond to a finite expectation value for the pion gradient in the direc- 
tion of the magnetic field, hence it will only be nonvanishing when chiral symmetry 
is broken. In that case, one has to extremize the on-shell action with respect to j 
ll70l l73l l74l . From the field theory perspective this means that, if j 7^ 0, the chiral 
condensate is rotating between a scalar and a pseudoscalar condensate when mov- 
ing along the z-direction, i.e., it forms a so-called chiral spiral l75l . Each full turn of 
the spiral raises the baryon number by one. Therefore, since j measures the rate of 
turns per unit length, it is related to the baryon density. Equivalently, one can regard 
j as a supercurrent, in analogy to superfluidity, where the phase of the condensate 
gives the superfluid velocity. 

Before continuing we put some restrictions on the gauge field and the embed- 
ding: in the joined configuration we assume that the fields are continuous at the 
junction point uq since for now we omit any point-like sources, hence 03(1*0) = 
and xfJuo) = 00 . In the restored phase, due to the presence of the horizon, we have 
to satisfy the regularity constraint flo( M r) = [67 1. 

Within our ansatz, the action for the D8-brane describing left-handed fermions 
becomes 



£l Juq/ut 



Vu 5 +b 2 u 2 ^ U lfxf + J- a' 2 + a' 2 h d - 
-y(a 3 a -aoa 3 ) 



(50) 



2 Here we keep the notation of Refs. [43 44) . Note that in the NJL section /i is the dimensionful 
quark chemical potential. 
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where the lower bound of the integration has to be chosen according to the phase 
under consideration and where 



T H R 5 Q 4 N C R 2 
67z 2 (2na' 



-= T^h- 7W - (51) 



Here we have modified the original action S and denoted the new action by S'. The 
reason is that proceeding with S results in an inconsistency: the conserved currents 
sourced by the boundary values of the gauge field turn out to be different from 
those currents calculated using the thermodynamic relations. In J73| this issue was 
related to the gauge variance of the CS action. The solution to this problem is to 
supplement the CS action with boundary terms residing at the holographic as well 
as at the spatial boundaries. After integration by parts this modification amounts to 
simply multiplying the CS contribution with a factor 3/2. 
The integrated equations of motion are 



Vu 5 + b 2 u 2 a' 



y/u 3 fx%-t 


- h -f-a' 2 + a' 2 h d 








- h f-a' 2 + a' 2 h d 




- b 2 u 2 fifix! 4 




-'f-a' 2 + a' 2 h d 



3ba^+c, (52) 



= 3ba + d, (53) 



(54) 



The left-hand side of Eq. ( 1521 is the magnitude of the (bulk) electric field corre- 
sponding to the gradient of ao in a curved background on one D8-brane pointing 
towards larger values of u. When we move past the point uq in the joined D-brane 
configuration the direction of the electric field is flipped since we assume that «o 
is P-even. Therefore, since 03 (uq) = 0, the integration constant c corresponds to a 
point-like source at uq. For now, we do not include any point-like baryons and thus 
set c = in the broken phase. In the restored phase, on the other hand, c ^ 0, hence 
the horizon provides a charge that will be translated into the quark density at the 
boundary. (In the restored phase, xfAu) = and thus k = 0, i.e., only two nontrivial 
equations remain.) Furthermore, if the magnetic field is nonzero there is an addi- 
tional contribution to the quark density from the gradient of 03, which in general is 
distributed over the whole D8-brane world volume. Equation ( 1531 evaluated at uj 
enforces us to set d = in the restored phase in order to maintain consistency since 
hd{ur) =a (u T ) =0. 

The nonvanishing components of the current densities sourced by the asymptotic 
gauge field components are given by 



2na'.yV f3b . 

y; + c), (55) 
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where we have used the equations of motion. The first line relates the baryon density 
with the magnetic chiral spiral and the point-like charges in the bulk. The second line 
is the axial current which we have already encountered in section [2~.2.2| Because we 
have to extremize the thermodynamic potential with respect to j, i.e., with respect 
to 03(0°), we can immediately conclude that in the broken phase the axial current 
has to vanish, hence d = — 3/2 bji. In the chirally symmetric phase the axial current 
at any temperature - reinstating dimensionful quantities - reads 

fI = ^Bn q . (57) 

This result differs from the corresponding expression (f28b obtained in the NJL 
model by a factor 2, which is related to the modification of the CS term in the 
action in order to obtain a consistent thermodynamic description of the currents. 
For a thorough discussion of the effect of this modification on the chiral anomaly 
seeRef. (761 . 



3.3 Semianalytic solution to the equations of motion 

In general, from this point on one has to rely on numerical methods. However, using 
f{u) ~ 1 we can go a little further. This approximation is valid either in the deconfin- 
ing background if T = or in the decompactified limit of the confined background. 
Moreover, as will be justified a posteriori, for L <C and in the chirally broken 
phase we have uq mkk (confined) or uq ^> itj for sufficiently small T (decon- 
fined). We will later work in the deconfined geometry and apply this approximation 
for the chirally broken phase at any T (i.e., our approximation becomes less accurate 
for large T). If chiral symmetry is restored this is of course not allowed, since the 
D8-branes extend from the holographic boundary down to the horizon at uj- Hence, 
when computing the phase diagram we will compare the grand canonical potential 
of the broken phase using the f(u) ~ 1 approximation with the full numerical re- 
sult obtained for the restored phase. Note that in the special case b = or jj. = 
the temperature can be easily introduced in the symmetric phase since the "black- 
ening" function f(n) does not appear explicitly in the equations of motion. There 
temperature enters only in the lower bound uj of the integrals over the holographic 
coordinate. 

With f(n) ~ 1, we can simplify Eqs. d52l ) and d53l considerably. We define the 
new coordinate field y(u) via the differential equation 



v/m 8 + u 5 b 2 -k 2 + (3b) 2 u 3 [{d y a ) 2 - (d y a 3 ) 2 } ' 
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sinhv M 


3^/icothVoo 









Table 1 The integration constants d, c\,C2,c,k and the supercurrent j for the chirally broken and 
restored phases. 



for which we have the freedom to choose y(wo) = or y(0) = in the broken and 
symmetric phase respectively. Its value at the holographic boundary will be denoted 
by Voo in the following. In the joined D8-brane configuration the boundary condition 
xfJuo) — > °° implies that y'(«o) — > °°. After algebraically rearranging Eqs. d52l - 
(T54-b such that all derivatives with respect to u are placed on the left-hand side, the 
equations of motion for the gauge fields as a function of the new coordinate y are 



dydo 



c 

3b' 



d y a3 = a 



3b 



(59) 



for which we can easily find the solutions 

d 



a 



■■ ci coshy + C2 sinhy - 



3b 



A3 = cj sinhy + C2 coshy - 



3b 



(60) 



This allows us to write the grand canonical potential, i.e., the on-shell action, as 



uq/ut y 



3b , kL 



^y- (4 - cf) - 1 (M - ci) + f (J - c 2 ) 



(61) 



This expression is divergent. In order to obtain finite expressions we renormalize 
the grand canonical potential by the chirally symmetric vacuum contribution 



Q(H = T = 0) = jY ( du\/u 5 +b 2 u 



(62) 



The integration constants found by imposing the boundary conditions discussed be- 
low Eq. d48l > and the supercurrent j = fl3(°°) are summarized in table Q] 



3.4 Broken chiral symmetry 

Inserting the supercurrent j and the constant c from table[T]into Eq. ( |55l l yields the 
quark number density 
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The only equations that remain and in general have to be solved numerically for the 
variables uq and y„ are 




y ( 3 ^ 
2 y -u 1 - o ^2coshy„ 



3b 



1W 



(64) 



where we have abbreviated 



(«); 




3b/i 



\2coshyo, 



b 2 4 



( 



\2coshyo. 



(65) 



Note that the explicit dependence on the asymptotic separation t can be eliminated 
by the rescaling u — > £ 2 u, jii — > l 2 \i, b — > £ 3 b and Q — > fQ. Therefore, in all plots 
shown below, the axes are measured in appropriate units of the D8-brane separation. 

Before coming to the full numerical results, let us first discuss the two limits of 
small and large magnetic fields b. For a detailed derivation of the approximations 
consult appendix D in Ref . B31 . 

For small magnetic fields, Voo and thus the supercurrent j rise linearly with b, and 
therefore the lowest order contribution to the quark number density induced by the 
chiral spiral is quadratic in b. The location of the tip of the connected flavor branes 
is uq ~ mI ' + Tji {pL)b 2 with the value of uq at b = 0, 



,(°) 



: 0.5249 . 



(66) 



Interestingly, the ji -dependent coefficient Tji possesses a zero at jtt ~ 0.2905/£ 2 , 
above which it becomes negative. This shows that the constituent quark mass (which 
is given by mo) can decrease with the magnetic field for sufficiently large chemical 
potentials. This behavior can be traced back to the incorporation of the chiral spiral. 

The grand canonical potential (renormalized by the vacuum contribution d62l )) is 
approximated for small b by 



ft, 



y(4 0) ) 7/2 



(67) 



where 



. 16/ 



K 

COt 

16 



.16) 



(68) 



(As explained in l43l , there exists a second solution in the region of small b, where 
wo is small and is large, which is separated from the solution discussed here by 
a first order phase transition. However, this first-order phase transition occurs in a 
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region of large /i where the chirally restored phase is preferred. Therefore we will 
not discuss this second solution here.) 

At asymptotically large magnetic field, diverges faster than linearly, thus j ~ 
jtt/2, while uq saturates at the value 



4y/7tr 



10 ' 



-i 2 



1.2317 £ 



-2 



(69) 



We see that > i.e., for any ji the constituent quark mass at asymptotically 
large b is larger than that at b = 0. This can be interpreted as magnetic catalysis and 
is similar to the NJL model. However, as we have shown in the left panel of Fig. 
|U in the NJL model the constituent quark mass does not saturate for asymptotically 
large magnetic fields. 

Plugging these results into £2 and n q yields 



On 



-jYb 



To 



2Tl 



Nc 

87T 



(70) 



Remarkably, all model parameters have dropped out of the quark number density, 
which thus is solely expressed in terms of the dimensionful quantities B and \i q . 



3.5 Symmetric phase 

The following analytical expressions are all valid in the zero-temperature limit. Only 
in the plots at the end of this subsection we include numerical finite-temperature 
results. Now only one equation remains to be solved numerically for y^, 

3bu 3 / 2 

=du. (71) 



For b > 0, this equation has in general three solutions: y^ = °°, which is always a 
solution, and two finite solutions, the larger of which turns out to be unstable. At 
sufficiently large values of b for a given pi only the divergent solution survives. For 
the quark density we find 

N c 

n q = ^B^cothjoo. (72) 

Let us first take the limit where b is small. In this case, yoo is linear in b, and we 
obtain for the (dimensionful) quark number density 
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VN C M KK 5/2 



5/2 



(73) 



The unusual exponent 5/2 of ji q can only occur due to the presence of the dimen- 
sionful model parameter Mkk (due to the extra dimension in the model), which 
provides the missing mass dimensions. 

The grand canonical potential becomes for small b 



On 




-,5/2 



(74) 



with 



'73 



5/2 



+ ■ 



r(|,)r(i) 



-M\ r 



. 10 ) 



(75) 



Taking the limit b — > °° allows only the solution >'oo = °°, as mentioned before. 
However, note that this is also a valid solution at finite b, hence the following results 
carry over to any value of b as long as this particular phase is considered. Interest- 
ingly, the density in this case is 



2^' 



(76) 



which takes precisely the form of the density of gapless free fermions in the lowest 
Landau level. Therefore, we may speak of a LLL-like phase in the Sakai-Sugimoto 
model, although there are, because of the strong-coupling nature, no quasiparticles 
and thus no Landau levels in the actual sense. The grand canonical potential is 



3bpL 2 



(77) 



Using (T74b together with d77b we can derive the critical magnetic field of the first- 
order transition within the chirally restored phase to the LLL-like phase as a func- 
tion of the chemical potential, 



b c ~ 0.095^ 



3/2 



(78) 



In the left panel of Fig. [H] we plot the quark number density for different tempera- 
tures. As a comparison, we also plot the corresponding density for (massless) free 
fermions in a magnetic field, obtained by taking the derivative with respect to the 
chemical potential of the thermodynamic potential (130b . 

In the case of free fermions, the higher Landau levels cause oscillations in the 
density at small magnetic field. These oscillations are absent in the "higher Lan- 
dau level phase" in the Sakai-Sugimoto model, given by the solution < °°. This 
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might be a consequence of the strong coupling, in which case we do not expect a 
sharp Fermi surface, even at T = 0. Furthermore, in the NJL model, the transitions 
between the phases with differently filled Landau levels, in particular also the tran- 
sition to the LLL phase, is second order, while in the Sakai-Sugimoto model it is 
first order. At finite temperature, the transitions become immediately smooth in the 
NJL model, while for given jj, it remains first order in the Sakai-Sugimoto model 
until a critical temperature is reached, which increases with increasing /i. Above 
this temperature only one minimizing solution for y„ exists for all b and given /i . 
As a result, the transition line in the b-pL plane has a critical endpoint for a given 
temperature, resulting in a critical line in the three-dimensional phase diagram, see 
Fig. [9] Another important difference is the location of the LLL-transition in the jj,- 
b diagram: the critical magnetic field at zero temperature is proportional to /i 3 / 2 , 
compared to /i 2 /2 for free fermions. Again this is due to the occurrence of a/Mkk- 




b 2B/^, 2 



Fig. 8 Quark number density as a function of the background magnetic field for a given chemical 
potential at various (dimensionless) temperatures f = TR in the Sakai-Sugimoto model (left) and 
the NJL model (right). 



3.6 Chiral phase transition 

First we discuss the critical temperature for chiral symmetry restoration at van- 
ishing chemical potential. In this case, in the restored phase the only tempera- 
ture dependence enters via the lower bound of the integrals over the holographic 
coordinate, uj = [Ant/3) , with t = RT. Therefore, one easily determines the 
renormalized grand canonical potential of the restored phase for the cases b = 0, 
Q mn = —2/ljVuj 2 , and b — > °°, Q Kn = — j¥ bvdj/1. Then, together with the corre- 
sponding expressions for the broken phase from Eqs. d67l l and d70l i we compute the 
critical temperatures 



t c [n=b = 0) =0.1355/^, 
t c {n =,Z>->°°) =0.1923/*. 



(79) 
(80) 
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(Remember that we have used the f(u) ~ 1 approximation for the broken phase 
which, strictly speaking, is only valid for very small temperatures.) We see that the 
Sakai-Sugimoto model reproduces the usual magnetic catalysis effect at zero chem- 
ical potential because the critical temperature at asymptotically large b is larger than 
that at vanishing b. This is supported by the numerical solution which shows that the 
critical temperature increases monotonically with the magnetic field. In contrast to 
the NJL model, the critical temperature saturates at the value given in equation d80t , 
because the value for mo, i.e., the holographic constituent quark mass, saturates. 

At zero temperature, we use Eqs. ( f67b and d70b for the broken phase and Eqs. 
d74b and d77l i for the restored phase to compute the critical chemical potentials 



This result already shows that inverse magnetic catalysis in the sense explained in 
Sec. l2.2.3l must be present in the Sakai-Sugimoto model. The full numerical solution 
of the surface of the chiral phase transition in the three dimensional T-pL-B space, 
including cuts through the surface at fixed t, jJ., and b, is shown in Fig. [9] 

In order to discuss the inverse magnetic catalysis, we have plotted the zero- 
temperature phase diagram separately in Fig. [10] This phase diagram shows in- 
triguing similarities with the corresponding NJL phase diagram in Fig. [6] inverse 
magnetic catalysis is present at small magnetic fields and is most pronounced when 
the restored phase has a LLL-like behavior. Even the manifestation of inverse mag- 
netic catalysis in the analytical approximations is qualitatively the same as in the 
field-theoretical model as we now show. For large magnetic fields, Eqs. (l70t and 
d74l i can be used to write the free energy difference between restored and broken 
phases as 



where we have identified Ruq/ (2na') with the constituent quark mass M Il66ll69l . 
This large-B expression for A £2 is remarkably similar to the weak-coupling expres- 
sion ( l38l in the NJL model. We can thus conclude, for the reasons explained below 
Eq. ( l38l ). that in the large-B regime the critical chemical potential must increase 
with B. This is confirmed by the chiral phase transition line of Fig. \TU\ Note the 
difference between the terms 6(jJ. — M) in the NJL expression and the last term 
in Eq. d83l l. Both terms come from a nonzero quark density which in our NJL cal- 
culation is only present if jx > M, while in our Sakai-Sugimoto calculation there is 
a topological quark density at nonzero B for all fx due to the chiral spiral. 

For small magnetic fields we may apply an approximation in the spirit of Eq. 
d40b . We compare the free energy of the broken phase for small magnetic fields d67l i 
with the free energy of the LLL phase ( ITTl i. The result can be written as 



H c (t=b = 0) = 0.4405 /f, 
y, c (t = 0,b^>-°°) = 0.4325 /£ 2 . 



(81) 
(82) 
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Fig. 9 Upper left panel: The surface of the chiral phase transition (blue) in the deconfined phase 
of the Sakai-Sugimoto model in the T-[l-B space. The small (green) surface shows the transition 
from the "higher LL" phase to the "LLL" phase, explained in Sec. 13.51 Upper right, lower left and 
lower right panels: two-dimensional cuts at various fixed temperatures, chemical potentials and 
magnetic fields, respectively, through the three-dimensional phase diagram. In the lower left plot, 
for instance, we see that the monotonically increasing critical temperature at fi = becomes a 
non-monotonic curve at finite [i and may even turn into two disconnected pieces, separating two 
chirally broken phases at small and large magnetic fields. 



where Mq °c is the constituent quark mass at B = 0. Again we recover the form of 
the NJL result (|40| |. The main conclusion is that the energy cost for condensation is 
linear in B, whereas the energy gain from condensation, i.e., the magnetic catalysis 
is only quadratic in B for small B. This allows for inverse magnetic catalysis. The 
dashed-dotted line in Fig. [10] is the approximate phase transition from Eq. ( f84T ). 
Comparison with the full numerical result shows that the approximation captures 
the physics of inverse magnetic catalysis where it is most pronounced and that the 
"hLL" phase counteracts inverse magnetic catalysis. 

In Fig. QT| we show the phase diagram including baryonic matter discussed in 
ll44ll . The main observations are that (i) baryonic matter prevents chiral symmetry 
restoration for small magnetic field for any value of n (as already found in Ref. l72l 
for B = 0) and that (if) for sufficiently large magnetic fields, baryons become disfa- 
vored, i.e., the chirally broken, mesonic, phase is directly superseded by the quark 
matter phase. Interestingly, in the presence of baryonic matter, inverse magnetic 
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Fig. 10 The chiral phase transition at zero temperature from the Sakai-Sugimoto model (ignoring 
baryonic matter). The chirally broken phase (white) is separated by a first-order phase transition 
(solid line) from the chirally restored phase (gray). The dashed-dotted line is the approximation 
from Eq. {84}. Translating the dimensionless quantities b and \i into physical units 1431 . one con- 
cludes that the magnetic field decreases the critical chemical potential from \l q ~ 400 MeV at 
\qB\ = down to fl c/ ~ 230MeV at \qB\ ~ 1 .0 x 10 19 G where the critical line turns around and the 
critical chemical potential starts to increase with \qB\. 
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Fig. 11 As Fig.[l0] but including baryonic matter (from Ref. 1441 ). The dashed line is the (second- 
order) onset of baryonic matter. The transition within the chirally restored phase between the "LLL" 
and "hLL" phases has disappeared because baryonic matter is preferred in this region of the phase 
diagram. 



catalysis becomes even more prominent in the phase diagram: now, the magnetic 
field restores chiral symmetry for any jj, > 0.25. 
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4 Discussion 

We have investigated equilibrium phases at finite temperature, chemical potential, 
and magnetic field for one massless flavor in the Nambu-Jona-Lasinio model and 
the Sakai-Sugimoto model. For small flavor brane separations, the Sakai-Sugimoto 
model is conjectured to be dual to a (non-local) NJL model. Indeed, we have found 
intriguing qualitative similarities between both models. 

There is an exact equality of the number density at zero temperature of the low- 
est Landau level in the restored phase of the NJL model and the large magnetic 
field phase with restored chiral symmetry in the Sakai-Sugimoto model. The higher 
Landau level phase in the NJL model, however, differs from the small magnetic 
field phase with restored chiral symmetry in the Sakai-Sugimoto model. For exam- 
ple, there occur no de Haas-van Alphen oscillations in the holographic model. One 
possible interpretation is that in the holographic model - dual to a strongly coupled 
gauge theory - there are no quasiparticles and no sharp Fermi surface. Furthermore, 
the axial current found on the field theory side is also reproduced in the holographic 
model. In the version of the model discussed here (73), the holographic current re- 
produces the field-theoretical current only up to a factor of 2. This discrepancy can 
be resolved by properly implementing the axial anomaly ll76l . however for the price 
of losing a consistent thermodynamic description. 

Also the phase diagrams in both models share the same qualitative features. The 
main differences are the order of the phase transitions (first and second order in NJL 
vs. first order in Sakai-Sugimoto), the saturation of the critical temperature and the 
critical chemical potential at asymptotically large magnetic fields (which only oc- 
curs in Sakai-Sugimoto), and the absence of de Haas-van Alphen oscillations of the 
phase transition line in the Sakai-Sugimoto model. The main physical effect, first 
discussed in detail in the holographic context ll43l . is the nontrivial behavior of the 
chiral phase transition in a magnetic field at finite quark chemical potential. Some- 
what unexpectedly, at sufficiently large chemical potentials and small temperatures 
and not too large magnetic fields, the effect of inverse magnetic catalysis dominates. 
We have explained inverse magnetic catalysis in both models by a free energy argu- 
ment. This argument shows that, even if the magnetic field increases the constituent 
quark mass (due to the usual magnetic catalysis) and thus increases the conden- 
sation energy, it also increases the energy cost for forming a chiral condensate. In 
particular, in the LLL, where the effect is most pronounced, the cost for overcoming 
the separation of fermions and antifermions due to the chemical potential increases 
linearly in B, while the constituent quark mass rises quadratically. It is interesting 
that at asymptotically large magnetic fields the free energy difference in the Sakai- 
Sugimoto model resembles the corresponding expression in the weak-coupling limit 
of the NJL model. In this regime magnetic catalysis is dominant in both models, 
and the situation is analogous to weak-coupling superconductivity with mismatched 
Fermi surfaces. 

By fitting the parameters of the holographic model with the help of the critical 
temperature at fi = B = from QCD lattice calculations |7J[8) and the (not very well 
known) critical chemical potential at T = B = from model calculations 117711781 . 
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we find that inverse magnetic catalysis persists up to B ~ 1.0 x 10 G, where the 
critical chemical potential has decreased from 400 MeV to about 230 MeV. It is not 
clear whether the magnetic field inside compact stars is large enough to have any 
effect on the chiral phase transition. Our results show, however, that if it is large 
enough then only inverse magnetic catalysis will play a role, i.e., the transition from 
hadronic to quark matter occurs at smaller densities than naively expected from the 
B = case. 

We have included an anisotropic chiral condensate in the Sakai-Sugimoto model, 
but not in the NJL model. For comparison, it is easy to show that in the holographic 
calculation the assumption of an isotropic chiral condensate does not change the 
qualitative features of the phase diagram. One finds that the effects of inverse mag- 
netic catalysis are rather enhanced. On the other hand, including an anisotropic chi- 
ral condensate in the NJL model changes the phase diagram drastically (79) ■ Most 
notably, there exists a phase with anisotropic chiral condensate even at B = 0; in the 
Sakai-Sugimoto model, B ^ is necessary for having such a phase. Moreover, this 
phase inevitably has a finite quark density. In order to realize this in the holographic 
model at B = one needs solitonic baryon sources which are related to Skyrmions 
and thus rather different from "baryons" in the NJL model which consist of dislo- 
cated quarks. We have briefly discussed the effect of such baryonic matter in the 
Sakai-Sugimoto model, based on Ref. 1441 . One of the most important changes is 
the non-existence of a chiral symmetry restoration at B = for any value of the 
chemical potential. 

Another phenomenon that was not included in our discussion is the so-called chi- 
ral shift (80 81 ], a chiral asymmetry in the Fermi surfaces of right- and left-handed 
charged fermions induced by a magnetic field. It would be interesting to discuss its 
effect on the chiral phase transition and thus on inverse magnetic catalysis. How- 
ever, the chiral shift is related to the Fock exchange terms, which are suppressed at 
large N c . Therefore, this effect is difficult to study in a holographic model where 
N c — > °° is necessary for the validity of the supergravity approximation. 

Acknowledgements This work has been supported by the Austrian science foundation FWF un- 
der project no. P22114-N16. 
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